Josephson super-current in graphene-superconductor junction 
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Within the tunneling Hamiltonian formulation for the eight-component spinors, the Josephson critical super- 
current has been calculated in a planar superconductor-normal graphene-superconductor junction. Coupling 
between superconductor regions and graphene is taken into account by a tunneling Hamiltonian which contains 
two types of tunneling, intra-valley and inter-valley tunneling. Within the present tunneling approach, we find 
that the contributions of two kinds of tunneling to the critical super-current, are completely separable. Therefore, 
it is possible to consider the effect of the inter-valley tunnelings in the critical super-current. The incorporation 
of these type of processes into the tunneling Hamiltonian, exposes a special feature of the graphene Josephson 
junctions. The effect of inter- valley tunneling appears in the length dependence plot of critical current in the form 
of oscillations. We also present the results for temperature dependence of critical super-current and compare 
with experimental results and other theoretical calculations. 

PACS numbers: 74.50.+r, 72.80.Vp 



INTRODUCTION 

In the past few years, after the first experimental synthesis 
by Geim et.al yj] the community has witnessed the growing 
interest in mono-layer carbon atoms arranged in a honeycomb 
lattice, called graphene. The peculiarity of graphene elec- 
tronic structure is the primary reason which has given rise to 
such enthusiasm among scientists. Linear dispersion relation 
near the discrete Fermi points and chiral nature of carriers in 
graphene are responsible for the most intriguing phenomena 
that have been reported 121 13|] . 

Among the attractive properties of graphene, from both 
basic and applied point of view, are the superconduc- 
tor proximity effect which has been studied experimen- 
tally and theoretically. Before any experimental observa- 
tion, Beenakker [|4J] predicted in a seminal work that at the 
graphene-superconductor interface, not only the conventional 
Andreev retro-reflection takes place, but also it can exhibit 
specular Andreev reflection owing to presence of two dif- 
ferent valleys in the graphene band structure. Specular An- 
dreev reflection takes place when an electron in the conduc- 
tion band converts into a hole in the valance band. At the 
low doping level, specular reflections are dominated and at 
Dirac point all reflections are of this type [Sfl. By solving 
the Dirac-Bogoliubov-de Gennes (DBdGD) equations for an 
"ideal" Normal-Superconductor (NS) interface, in the short 
junction limit, in which coherence length ^ = hvp/^o {vf 
is the Fermi velocity and Aq the superconductivity order pa- 
rameter), is much larger than junction length, L, and with ne- 
glecting inter-valley scatterings, Titov and Beenakker ||6[] and 
simultaneously Moghaddam and Zareyan 17|] calculated the 
Andreev bound states, with which they obtained Josephson 
super-current. Their results predicted the existence of a finite 
cuiTent even at the Dirac point. In these works, to solve the 
DBdG equations, the rigid boundary condition was assumed, 
according to which superconductivity gap has a fixed and fi- 
nite value in the superconductor regions and it is zero in the 



normal region. 

First experimental investigation of superconductor- 
graphene-superconductor (SGS) Josephson junction jsl] 
showed that Josephson current does flow through these junc- 
tions. This current depends strongly on the position of Fermi 
energy and as it had been predicted, there was a non-zero 
super-current at the neutral Dirac point. Afterwards, several 
experimental studies have been conducted i9t4l3ll. Recently, 
in addition to critical super-current, direct measurement of 
the current-phase relation (CPR) has been performed with the 
interferometry technique il4l llSll . The results confirm the 
predicted non-sinusoidal CPR curves 161] and show that the 
deviation from the sinusoidal behavior increases linearly with 
critical current. 

At the same time, more theoretical efforts were trig- 
gered I15I422II . Self-consistent solution of the tight binding 
DBdG equation, at zero temperature JloL llTIl and later gen- 
eralization to finite temperature llSll . have provided results 
on the dependence of the Josephson super-current on length, 
temperature, doping level, phase difference and pairing sym- 
metries of the superconductivity order parameter Another 
method based on Cooper pair propagation over long distances 
for junctions with L ^ W also reveals such dependencies in 
the graphene Josephson junctions il9ll20ll . Shifting the Fermi 
energy away from zero point results in an enhancement in the 
critical super-current due to finite density of states at the Fermi 
level for /i 7^ 0. Besides the enhancement, one can see oscilla- 
tions in the diagrams of the critical super-current as a function 
of length JilOl- 

In this paper we will use the perturbative Green's function 
method in framework of path integral and tunneling Hamil- 
tonian between superconductor and normal graphene regions 
to calculate the critical Josephson super-current. An s-wave 
superconductivity pairing will be assumed in the supercon- 
ductor areas of the graphene honeycomb lattice. In contrast to 
previous works, we don't neglect the inter- valley processes, 
as such processes may become important in the presence 



of edges, or external potentials with sharp variations 123 1. 
Among the other circumstances where inter-valley scattering 
becomes imp ortant, one can mention the presence of Kekule 
textures 12411 which may arise from creation of instantons of 
opposite sign in two valleys 12511 . We will incorporate both 
types of inter-valley and intra-valley tunneling processes into 
the Hamiltonian. We will discuss the temperature and length 
dependence of the critical super-current. The advantage of the 
present tunneling approach is that, within this approach, con- 
tributions of inter-valley and intra-valley processes to the total 
Josephson super-current can be totally separated. Therefore, 
the effect of the inter- valley tunneling on the critical current 
can be isolated. The present formulation reveals a peculiar 
feature of the graphene Josephson junctions which is the re- 
sult of incorporating the inter-valley tunnelings. In the curve 
of critical super-current as a function of junction length, the 
effect of inter-valley tunneling appears in the form of oscilla- 
tions which are drastically different from oscillations due to 
non-zero chemical potential. 

The paper is organized as follows. First, we summarize 



the perturbative Green's function method 12611 and rephrase it 
for the 8-component spinors needed in the present paper. We 
then apply the formulation to the problem of SGS junctions. 
We close the paper by a discussion on the results. 



METHOD AND MODEL 

In a Josephson junction between two superconductors, a 
dissipation-less super-current flows from one superconductor 
to the other to provide a state with minimum energy 12711 . The 
critical super-current is given by. 



T 2e^i^ T ■ I \ 



(1) 



where F is the free energy of the system, Ic indicates the crit- 
ical super-current and ip = (pn — ipL is phase difference be- 
tween two superconductors. In the framework of the path in- 
tegral, partition function is given by. 
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in which S is the effective action of system that is a function 
of Grassmann variables ^ and ^. In the imaginary time for- 
maUsm (r = it), the action will be given by, 

5 = V / dr[*(/c,r)9^*(/c,r) + H(*(fc,T),*(fc,T))]. 

(3) 
Here, H is the full Hamiltonian of the system that takes the 
form 



H = Hl + Hr + Hn + Ht. 



(4) 



Hl, Hfj, Hpf and Ht are respectively the Hamiltonian of the 
left and right superconductors, normal region and tunneling 



between superconductors and normal regions. The form of 
these Hamiltonians can be quite general. So that the for- 
mula ( fT2l i holds under very general circumstances. The par- 
ticular form of these Hamiltonians suitable for our own prob- 
lem are given in Eqs. ( fT3b . ( fT4l i. ( fTSl ). If we express the 
matrix form of the above Hamiltonian as sum of two parts, 
T as tunneling term and g^^ for the rest, then in the basis 
#(fc) = (>l'(fci,),*l'(fcAr),'l'(fcfl,)), it can be written in the 
compact matrix form as, 

H = J2'^\k)[g^'+fMk). (5) 

{k} 

where ^q"^ and T are given by the following matrix forms: 
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Here, g'^^, g^^ and g^^ are respectively the inverse propaga- 
tors for the left, normal and right areas excluding the tunneling 
parts. 

After a Fourier transform, the effective action becomes. 



S=Y1 *(fc,iw„)[-iw„I + 5o^+T]*(fc,ia;„ 



(8) 



where hujn = 7r/sBT(2n + 1) are the Matsubara frequencies. 
With definition. 
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and from Eq. (|2|i we have, 

(10) 
This is a Gaussian integral, which can be performed to lead to 
the following equation for the free energy, 

1, 

13' 
1, 



F=-^Trln[-Goi+f], 
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Tr[ln(-Goi) + ln(l-Gor)]. (11) 



Here, Tr means a summation over all diagonal matrix ele- 
ments and an integration over all momenta. Having expanded 
this expression in terms of the tunneling amplitude, it can 
be readily shown that the fourth-order term is the first non- 
zero and leading term that contributes to the Josephson super- 
current. Using Eq. ([T]i and the identity In dot A = Tr In A, the 
final expression for super-current is obtained as, 

— -Tt[GlTnlGnT^rGbTnb.GnT^i^]. 
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FIG. 1: The schematic geometry of graphene Josephson junction. L 
is the separation of two superconductor (Sc) regions and W is the 
junction width. 



Now let us specialize the general formula ( fT2] ) to the case 
special case of SGS junctions of this paper The geometry 
of the graphene Josephson junction that we consider in this 
paper is depicted in Fig. [T] In the normal graphene region, 
the low energy electrons are governed by the relativistic Dirac 
Hamiltonian, Hn = vp{kN ■ a) (E) tq, where vp — lO^m/s is 
the Fermi velocity in graphene and a and t are Pauli matrices. 
a operates on the sub-lattice degrees of freedom and tq acts 
on the valley degree of freedom. Dirac fermions in graphene 
live in two different valleys, K and K, on the opposite corners 
of the Brillouin zone. 

Pristine graphene cannot be superconductor in ordinary 
conditions. Experimentally, to prepare an SGS junction, 
superconducting electrodes can be deposited on top of a 
graphene sheet, so that due to this proximity of supercon- 
ductor electrodes to graphene layer, the superconductivity or- 
der parameter can be assumed to be induced on the left and 
right graphene regions. Here, an s-wave superconductivity 
will be assumed in the left and right superconducting areas. In 
our model, the assumption is that the superconducting pairing 
takes place between two time-reversal electrons on the same 
sub-lattice but on different valleys in the Brillouin zone 12811 . 
Thus the mean field Hamiltonian for the superconducting part 
will be. 
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(13) 
where a and 6 are fermion operators on sub-lattices A and B 
pertaining to valley K, while a and 6 denote the corresponding 
operators on valley K. Aq is the magnitude of superconduct- 
ing order parameter and (pa denotes phase of each a = R,L 
superconducting leads. 

The electron tunneling between superconductors and nor- 
mal graphene region can be considered in several different 
ways. In our case, as it is depicted in Fig. [T] superconduc- 
tor and normal graphene are assumed to be connected to each 
other through their zigzag edge. Tunneling Hamiltonian con- 
tains left and right tunneling between superconductor and nor- 
mal graphene. As can be seen in Fig.[Tl for right side, electron 
tunneling takes place between sub-lattice A in normal region 
and sub-lattice B in superconductor region while at the left 



interface, electrons tunnel between sub-lattice B in normal re- 
gion and sub-lattice A in superconductor area. Therefore, the 
tunneling Hamiltonian in tight binding approximation can be 
written as. 
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(14) 



where jt is the tunneling amplitude and A and B generally 
stand for the Dirac fermion operators on the A and B sub- 
lattices and summation extends over all atomic sites along the 
interfaces. By writing the Dirac fermion operators as the sum 
of two valley operators, tunneling Hamiltonian for the right 
side will be given by. 
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and likewise for the left interface. As can be seen, tunneling 
Hamiltonian contains two types of tunneling, intra- valley (first 
and last terms) and inter-valley (second and third terms) tun- 
nelings. Inter-valley tunneling involves a "momentum trans- 
fer" Q, that connects the two independent valleys in the Bril- 
louin zone. Now, let us define an appropriate basis to rewrite 
all parts of Hamiltonian. A suitable representation for Fermi 
operators is the following: 
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In this basis, Gq^ and GqJ (a = R, L) are block diagonal 
and one can easily obtain the required matrices by just invert- 
ing the blocks. The Green's function of the normal region is 
written as. 
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(17) 
It is important to mention that in the normal area two electrons 
propagate between two superconductors as a Cooper pair with 
different spin. Thus, kN-^ and k^^ are independent degrees 
of freedom in this region. For the superconductor Green's 
function we have. 



GqI = *^rJi 
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Finally, for tunneling part of the Hamiltonian we have. 
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where we use these notations: 
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(21) 
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Now, we have all of matrices needed for calculation of the free 
energy. From Eq. (fT2] i and by using these matrices, we obtain: 
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F ^- 



16A^j^{hvF)^ [(1 + cos9)kikl + (-1 + cos9)k^ykl - siii0(44 + 4^)] cos((^ + (^ - fct)i) 
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where A: is a collective index denoting all possible momenta, 
i.e. k — {kL,kR,klf,k]^}. In Eq. ( |23] ), we have used the 
following notations. 
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The summation in Eq. (l23T l is over all Matsubara frequencies 
and all momenta around the pertinent valleys. In the limit of 
wide junction limit, i.e. W ^ L, the details of the interface 
becomes UTelevant and we can replace the summation with an 

I 



integration around the corresponding valley, up-to an energy 
cutoff that the linear dispersion holds. 

Due to symmetry of the integration region, even terms of in- 
tegrand can contribute to the free energy. With definition ([T]i 
for critical super-current and from Eq. ( |23] l we obtain (see Ap- 
pendix), 



, l28ekBTA^{Thl.SHl + cosQ ■ L) ^ ^ ^^. : 

Ic = -^ 2^ [f{lUJn)g{lUJn)\ 

(25) 
where S — WL is the area of junction and functions / and g 
are given by. 
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(27) 
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Ec is the cutoff energy and L indicates the vector connecting 
two superconducting leads. 

Carefully tracing the above derivation, shows that in 
Eq. ( |25] ), all effects of the inter-valley tunnelling processes 
appears in a term proportional to cos((5 • E). The elimi- 
nation of the term proportional to sin(Q • E) can be physi- 
cally understood: In view of symmetry of the problem un- 
der y ^ —y transformation, the Hamiltonian is symmetric 



under fc„ 



-ky, which imphes that free energy must be 



symmetric with respect to Q -^ —Q transformation. Hence 
only the term proportional to cos(Q • E) will survive in the 
free energy. Therefore, an oscillatory behavior in the crit- 
ical super-current is predicted which is solely due to inter- 
valley processes. It can be seen in the final result by setting 



Q = 0, which amounts to ignoring the inter-valley processes, 
(1 + cos{Q.L)) -^ g-i, = 2. This is the simple valley degener- 
acy expected in all extensive quantities, when the inter-valley 
processes are ignored 12311 . 

The fact that contribution of inter- valley processes (both in 
tunneling and superconducting pairing) gives rise to an addi- 
tive modulation by a wave vector Q is a feature of the present 
tunneling formulation. This shows that within the present 
formulation the contributions of two types of tunneling to 
the critical super-current are naturally separable, and repro- 
duces to proper degeneracy factor g^, = 2 in the limit where 
inter-valley processes are ignored. The oscillatory factors of 
1 + cos(Q.r) arise in other contexts, which require a convolu- 
tion of Green's functions related to Dirac cones 12911. 
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FIG. 2: Plot of the critical super-current as a function of junction 
length for T — O.lTc. The solid line with sharp oscillations shows 
the total critical current and the dashed line presents the contribution 
of the intra- valley tunneling processes to the critical super-current. 



RESULTS AND DISCUSSION 
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FIG. 3: The critical super-current as a function of junction length 
for different temperatures. For easier comparison, in this plot, the 
inter-valley tunneling processes have not been taken into account. 
Decreasing the critical super-current by increasing junction length 
and temperature, is clearly seen in these curves. 



To calculate the critical super-current from Eq.jZSll we need 
to have the temperature dependence of the magnitude of su- 
perconductivity order parameter, AfT). For s-wave supercon- 
ductors the pair potential satisfies OOll . 



In 
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where Ko is the modified Bessel function and Ao(T = 0) 
is the magnitude of the superconductivity gap at zero tem- 
perature. Analytical calculation of Eq. dZSl l is not feasible, 
thus we will calculate it numerically. Electron dispersion re- 
lation in graphene remains linear below a momentum cutoff 
that coiTesponds to Ec = hvpkc ~ leV. Careful scrutiny of 
the Eq. ( |27] | reveals an lj~^ dependence of the summand, so 
that it would rapidly decay at larger Matsubara frequencies. 
Therefore, it is adequate to sum over only first few frequen- 
cies. We are now ready to present the results of calculations 
of the critical current using Eq. dZSl l. As indicated before, we 
are in the wide junction limit {W ^ L), but we have no re- 
striction on the distance between the superconductors, i.e we 
can calculate the critical current for all range of the distances 
(L < ^ , i ~ ^ and L > ^). In all of our calculations we set 
W = 10/im > C 

Figure |2] shows the length dependence of the critical cur- 
rent. The temperature at which the calculation has been per- 
formed is O.lTc- In this figure, the curve with sharp oscilla- 
tions shows the total critical current, while the smooth digram 
gives the contribution of the intra-valley tunneling processes 
to the critical current. As it is obviously expected, the cur- 
rent decreases by increasing the distance between two super- 
conductors. Although the decreasing trend of two curves is 
similar to each other, but the oscillatory behavior of the total 
critical current discriminates between two types of tunneling. 
In fact, the distinguishing feature of the inter-valley electron 



tunneling is the oscillation with a characteristic wave-vector 
IQI in the total current. This suggests that if one measures the 
critical current versus the junction length, an oscillatory be- 
havior is expected which can be attributed to the inter-valley 
tunnelings. 

The presence of a second valley connected to the first 
one by time reversal symmetry, is the reason for emergence 
of a novel specular Andreev reflection 14j, |5[] . The present 
tunneling formulation reveals that existence of two valleys 
leads to another peculiar feature in Josephson current through 
graphene, namely oscillatory behavior of the Josephson cur- 
rent as a function of the junction length, L. In fact, our results 
in Fig.|2]shows that for some distance between two supercon- 
ducting electrodes, the critical super-current will be strongly 
suppressed due to inter-valley tunnelings between normal and 
superconducting regions. This can be physically interpreted 
as the destructive interference between the super-current aris- 
ing from Andreev bound states |]7[], and those arising from An- 
dreev modes [4J] (made possible by specular reflections). Note 
that, in all our calculations, an undoped normal graphene re- 
gion has been assumed, so that, undoubtedly the doping 11191 - 
2111 is not responsible for the oscillatory behavior of the criti- 
cal super-current. 

Figure [3] presents the length dependence of critical super- 
current for different temperatures. In this plot, to show the 
general dependence of the critical current on the junction 
length more clearly, and also to facilitate comparison with 
other works which do not consider inter-valley tunneling, we 
have only retained the intra-valley tunnelings. It is clearly 
seen that the critical super-current suddenly decreases in all 
diagrams, when one approaches L k, ^ from below. Fur- 
thermore, regarding the temperature dependence, it drops very 
rapidly when temperature approaches the critical temperature. 
Our results are in a good agreement with the results of the 
self-consistent tight-binding method (Fig. 5a in Ref. lll6ll and 
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FIG. 4: Temperature dependence of the critical super-current for dif- 
ferent junction length (up) for L < ^ and (down) L > ^. Notice that 
in these diagrams only the contribution of the intra-valley tunnelings 
to the critical current, have been taken into account. 



Fig. llainRef. |ll7|]) 

In Figure |4] we have plotted the contribution of the intra- 
valley processes in the critical super-current as a function of 
temperature. The results are displayed separately for distances 
below the coherence length, L < ^ (top panel) and distances 
above it, L > ^ (bottom panel}. The top panel of this fig- 
ure agrees with results of Ref. 12011 . When one compares the 
temperature dependence for junctions lengths less than the co- 
herence length, with those above the coherence length, a clear 
qualitative difference can be observed: For distances below 
the coherence length, there is a plateau in super-current for 
temperatures smaller than 0.5Tc and after that all diagrams 
reach zero near to T = Tc- But for distances above the co- 
herence length, after a much shorter plateau, an exponential 
decay takes place for temperatures near Tc- As one can see 
in the Fig.|4] the temperature at which rapid decrease sets in, 
decreases by increasing the junction length. 

Temperature dependence of the critical current has been 
measured in several research works. Our results are qualita- 
tively consistent with these experimental measurements. In 
particular, when a comparison between the observations in 
Ref. lull and Ref. lll3ll is made, the above-mentioned different 
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FIG. 5: Plot of the total critical super-current as a function of tem- 
perature for different values of the distance between superconducting 
electrodes. In contrast to the diagrams of the Fig. |4l in this diagram 
the contributions of all processes have been calculated. 



behaviors of temperature dependence of the critical cuiTent 
can be noticed. Therefore from our calculations we judge that 
in samples of Ref. lllll . the junction length is smaller than the 
coherence length, while for those in Ref. 11311 . the length of 
the junctions seems to be above the coherence length. Similar 
temperature dependencies have been obtained in the other the- 
oretical works. For example, in Fig. 2 of Ref. jlSll and Fig. 3 
of Ref. 120(1 , for L < ^ and L > ^, similar results have been 
obtained. Finally, in Fig.|5]the total critical super-current as a 
function of temperature has been plotted. In this diagram both 
inter-valley and intra-valley tunnelings have been taken into 
account. As can be seen in this figure, the variation of critical 
current with the junction length is different from diagrams of 
Fig.|4] which is a consequence of the oscillatory behavior the 
critical current as a function of the junction length. Otherwise, 
the general trends are similar to those in Fig.|4] 

In summary, we employed a method based on perturbative 
Green's function in the framework of path integral, to calcu- 
late the critical super-current in the graphene Josephson junc- 
tion. Our results presented the length and temperature depen- 
dence of the critical super-current. The role of the inter-valley 
tunnelings in the super-current has been investigated and it is 
shown that incorporating these kind of tunnelings, led to sharp 
oscillations in the behavior of critical current as a function of 
the junction length. Comparison of our results with the results 
that have been obtained by the other formalisms and also with 
experimental observations, suggests that the junctions realized 
in experiment are likely to be in weak tunneling regime which 
are consistent with perturbative treatment. 
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Details of calculations 

To obtain Eq. dZSl ) from Eq. ( l23b , we should expand the nu- 
merator of Eq.|23] For convenient, we replace k^L, k^L, k^L 
and k^L respectively with a;, x\ y and y'. There are two main 
terms: For the one proportional to sin cp, we have, 



sui(/? 



(1 + cos6')a;a;' + (— 1 + cos 6)yy' — smd{xy' + x'y) 

(29) 



since cos x — cosccsma; 



These terms are odd functions at least with respect to one of 
their variables while the integration region is even respect to 
all variables, hence they do not have contribution to the free 
energy. The other term is proportional to cos (p, 

costpi (1 + cos 9)xx' + (— 1 + cos 0)yy' — sin 8{xy' + x'y) 

[ cos a; cos a:' + sin x sin x'] > . (30) 

After careful consideration we find out that only the term pro- 
portional to X x' sin X sin x' has non-zero contribution to the 
free energy. Therefore, the remained integration correspond- 
ing to normal region variables is of this type: 



dy I dx 



Ea 



xsmx 



ygp^ {hwnY + x^ + y'^ 



(31) 



Integration of variable y can be performed analytically, which 

gives. 



X SVTLX 



2 dx 

Jo ^/{hLOn)^+X^ 



arctan - 



E'^. 



{hiuny 



(32) 
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